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Degeneracy in the characterization of non-transiting 
planets from transit timing variations 
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ABSTRACT 

The transit timing variation (TTV) method allows the detection of non-transiting 
planets through their gravitational perturbations. Since TTVs are strongly enhanced 
in systems close to mean- motion resonances (MMR), even a low mass planet can pro- 
duce an observable signal. This technique has thus been proposed to detect terrestrial 
planets. In this letter, we analyse TTV signals for systems in or close to MMR in order 
to illustrate the difficulties arising in the determination of planetary parameters. TTVs 
are computed numerically with an n-body integrator for a variety of systems close to 
MMR. The main features of these TTVs are also derived analytically. Systems deeply 
inside MMR do not produce particularly strong TTVs, while those close to MMR 
generate quasiperiodic TTVs characterised by a dominant long period term and a low 
amplitude remainder. If the remainder is too weak to be detected, then the signal is 
strongly degenerate and this prevents the determination of the planetary parameters. 
Even though an Earth mass planet can be detected by the TTV method if it is close 
to a MMR, it may not be possible to assert that this planet is actually an Earth mass 
planet. On the other hand, if the system is right in the center of a MMR, the high 
amplitude oscillation of the TTV signal vanishes and the detection of the perturber 
becomes as difficult as it is far from MMR. 
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1 INTRODUCTION 

A transiting planet belonging to a multiplanetary sys- 
tem may suffer gravitational interactions. The perbuta- 
tions in the orbit can then produce observable variations 
in the transit midtimes with respect to those of a purely 
Keplerian orbit. This is the basic idea behind the tran- 
sit timing variation method which has been proposed to 
detect non-transiting te rrestrial planets l|Agol et al.l 120051 ; 
iHolman fc Murray! l2005h . 

The TTV si gnal scales roughly linearly with the mass 
of the perturber (lAgol et al.ll2005l ; IHolman fc Murrav||2005l ; 
iNesvornv fc Morbidellil 120081 ). Thus, low-mass planets gen- 
erate weak TTVs except if they are in mean motion reso- 
nance o r near period commensurability with the t ransiting 
planet |Agol et al. I l2005l ; IHolman fc Murray! I2005D . In that 
case, the TTV signal contains a long term oscillation as- 
sociated to_a_^gj[eat i j2ieajiah 

system (|Lissauer et al.ll201ll ; [Ballard et alj|201ll h This long 
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period, like all the oscillations in the TTV signal, scales lin- 
early with the period of the transiting planet and is of the or- 
der of a year for a three-day transiting planet. Such systems, 
if they exist, should thus start to be detectable in surveys 
like CoRoT or Kepler. However , most of the systems close 
to MMR generate similar TTVs jFord fc Holmanll2007l 'l. and 
many p erturber parameters can reproduce an observed TTV 
signa l (TNesvornv fc MorbideTiH |200S| ; iMeschiari fc Laughlml 
2010; IVeras et al.ll201l| y Such a difficu lty has already bee n 
experienced with the Kepler 19 system (|Ballard et al.ll201ll ). 
Nevertheless, short-term perturbations, if detectable, should 
produce a "choppin g signal" that raises the degeneracy 
(|Holman et al.ll2010ri . 

Analytical expressions of the order of magnitude of the 
amplitude and the period of TTV signals for systems in 
or near first-orde r mean-motion resonances are provided by 
lAgol et al.l(|2005l h Here, we propose to focus on higher order 
resonances. In a first section, we compare the TTV signal 
produced by an Earth mass planet with typical detection 
thresholds to motivate the study of these resonances. In a 
second section, we present a few examples of TTV signals for 
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a variety of planetary parameters. In the subsequent section, 
we provide analytical approximations of the TTV signal for 
systems near or at MMR. We conclude in the last section. 



2 TTV AMPLITUDE PRODUCED BY AN 
EARTH MASS PLANET 

We consider the strength of TTVs produced by a terrestrial 
planet of mass m,2 = M$ on a Jupiter mass planet (mi = 
Mj) transiting a star with mass mo = M©. The orbit of the 
transiting planet is assumed to be initially circular with a 
period Pi = 3 days. Since second and higher order mean 
motion resonances only exist for non zero eccentricities, we 
set the eccentricity of the perturbing planet to 62 = 0.1. 
The TTVs are computed using an n-body integrator over 
N = 300 transits which corresponds to 2.5 years. The system 
is supposed to be coplanar, the initial mean longitudes Ai, 
A2 of the two planets and the longitude of the periapsis W2 
of the perturber are arbitrarily set to 270 deg which is the 
direction of the observer. 

The strength o f TTVs is measured in the same way as 
m lAeol et al l fijooih. using the standard deviation 
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where Pi and to are chosen to minimize a, and tj, j ^ 1, are 
the midtransit times of the N transits. This quantity a is dis- 
played in figure [1] for 1000 different initial conditions evenly 
spread over the period P2 of the perturber between 2Pi 
and lOPi. We recognize the well known peaks centered on 
period commensurabili ties corresponding to amplified grav- 
itational interactions l|Agol et all 120051 : I Holman fc Murray! 
2005). The figure shows also the range of typical detection 
thresholds varying from 10 sec to 1 min depending on the 
instrument and on the depth of the transit. It is notewor- 
thy that, for the system considered in this letter, terrestrial 
planets can hardly produce signals at the limit of detection 
and only in the vicinity of period commensurabilities with 
P2 /Pi ^ 3. Since the amplitude of TTV scales linearly with 
the period of the transiting planet, the signal produced on 
a wider orbit, with a period ten times larger for instance, 
would be easier to detect. However, the period of the TTV 
signal would be ten times longer too. It may exceed the 
duration of the surveys like C0R0T and Kepler. With this 
choice of initial conditions, most of the systems do not enter 
MMR because the value of the resonant angle is that of the 
hyperbolic point (see section l4~2j) . 



3 TTV SIGNAL NEAR A PERIOD 
COMMENSURABILITY 

To illustrate the degeneracy of the TTV signal, we consider 
8 different initial conditions, near different mean motion res- 
onances. The system is the same as in the previous section, 
apart from the mass of the perturber m,2 and its eccentricity 
e2 which are chosen so as to produce TTVs with the same 
amplitude (~ 2 min) and the same period (4 oscillations 
during the 2.5 years) for each of the different orbital periods 
P2. The initial conditions are gathered in Table Q] Although 
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Figure 1. Standard deviation of TTVs produced by an Earth 
mass planet on a Jovian planet transiting a Solar mass star every 
three days. The stripe shows the range of typical detectability 
thresholds. 



Table 1. Planetary parameters used to perform the simulations 
of the Figure [2] 



set a 


MMR 


m 2 (Af ffi ) 


e 2 b 


£ 




1 


1 : 2 


0.9 


0.087 


0.0132 


1.36 


2 


4 : 9 


24.5 


0.120 


0.0034 


1.15 


3 


3 : 7 


21.1 


0.100 


0.0044 


1.13 


4 


2 : 5 


8.6 


0.102 


0.0067 


1.02 


5 


3 : 8 


17.2 


0.160 


0.0045 


1.07 


6 


1 : 3 


49.7 


0.100 


0.0134 


1.25 


7 


3 : 10 


95.2 


0.194 


0.0045 


1.30 


8 


1 : 4 


394.0 


0.115 


0.0134 


1.31 



a The initial conditions are in order of increasing semi-major axis. 
b The eccentricity e2 is chosen to produce low values of the re- 
duced chi square \fx%.. 



this letter focuses on second and higher order period com- 
mensurabilities, we consider the 1:2 resonance to highlight 
the fact that even a 0.9 Earth mass planet can generate a 
TTV signal with the same amplitude and shape as a 4.14 
Saturnian mass planet near the 1:4 resonance. 

Figure [2] shows the TTV signals, computed for the 8 
initial conditions of table Q to which a gaussian noise with 
Cnoisc = 20 seconds has been added to simulate observa- 
tions. All the signals are very similar and easily detectable. 
But they are dominated by one single frequency and look 
like simple noisy sinusoids. It is thus quite difficult to differ- 
entiate them. For each signal, the highes t amplitude s i nusoi d 
is extracted using a frequency analysis ( La skar|[l99d . ll99o^ ) . 
and represented by a white curve in Fig. [2] Then, the re- 
duced chi-square of this fit is computed and shown in the 
last column of the table [T] In all cases, the fit by a pure sine 
leads to a reduced chi-square lower than or of the order of 
1.30. Naturally, if the two planets transit, the degeneracy 
disappears because the resonance is known. Here, we focus 
only on non-transiting perturbers. 

In the following we derive an analytical approximation 
of this highest amplitude sinusoidal term. 
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Figure 2. Example of TTV signal produced near period com- 
mensurabilities. A gaussian error of 20 seconds has been added 
to the signal to model the observation uncertainties. The white 
curve represents the highest amplitude term obtained by fre- 
quency analysis. The gray curve is the analytic approximation 
Eq. 121> . The number at the upper right corner of each subfigure 
corresponds to the set of initial conditions as given in table [T] 



4 ANALYTICAL APPROXIMATION OF THE 
HIGHEST AMPLITUDE TERM 

We consider a coplanar system in the vicinity of the p : p + q 
mean motion resonance. We assume that the orbit of the 
transiting planet is initially circular, while that of the per- 
turber is eccentric. We define Oj, Ai, e;, Wi as the semi-major 
axis, the mean longitude, the eccentricity and the longitude 
of pericenter of the ith planet, respectively. The conjugated 
variable to the mean longitude Xi is, up to the first order 
in planet masses, A, = rriiy/Gmoai with G being the grav- 
itational constant. The index i = 1 stands for the inner 
(transiting) planet, and the index i — 2 refers to the outer 
perturber. For each orbital parameter Ei, we note the un- 
perturbed part Ei, and the oscillatory perturbation SEi such 
that Ei — Ei + SEi. For the system analysed in this section, 
the terms contributing to the transit timing variations St of 
the inner planet are (e.g. iNesvornv fe Mor bidclli 200i|) 

mSt = d\i + 28k! sin(Ai) - 28hi cos(Ai) + 0(ei) . (2) 

In this expression, m = 2w/Pi, defined by Gmo = n\a\, 
is the mean motion of the transiting planet, ki = ei coswi 
and hi = e\ sin w\ are the usual non-singular eccentricity 
variables of the transiting planet. With the convention Ai = 
37r/2 at transit, Eq. ([2]) becomes 

mSt = 8Xi - 28ki + 0{e x ) . (3) 



The Hamiltonian H governing the evolution of a copla- 
nar two-body system can be decomposed into a Keplerian 
part 

Gmomi Gm§m 2 



2ai 2<J2 
and an interaction Hi given by 

Gra\m,2 



Hi = 



a-2 



(4) 



(5) 



x cos (ji Ai + J2A2 + + J4ZU2) , 

where a = 01/02. C](q) with j = (ji, J2, J3, J4) and 1 = 
(^1,^2) are coefficients that can be computed either in terms 
of Laplace coefficients b^ 2 (a) (e.g. [Ellis fc Murravl [20001 ) . 
or through a series in power of a (e.g. Kaula 1962). 

Here, we make the assumption that the evolution of a 
system in the vicinity of the p : p + q MMR is well recov- 
ered when only the resonant terms with the lowest power 
in eccentricity are kept in E q. j5l). This is verified if the 
resonances are well separated (|Wisdoml ["19801. One gets 



(6) 



Hi ~ —e\Hi cos (pAi — (p + g)A 2 + qzu 2 ) 

— eiel~ 1 H 2 cos (pAi — (p + g)A 2 + m\ + (q—l)zu' 2 ) 
where Hi = 2Gmim 2 Ci{a) / '02 and 

Ci(a) = C{° : %_ qM (a) , C 2 (a) = Cfr^^a) . (7) 

The factor 2 in the definition of Hi comes from the symmetry 
C](q) = C_j(a). In JBJ, we consider only the linear terms in 
ei as the higher orders do not contribute to the equations of 
motion for ei = 0. Since the eccentricity e 2 does not enter 
in the expression of St @, we make the hypothesis that e 2 
and its longitude of pericenter vj 2 are dominated by their 
constant (initial) value for the period of time considered 
in this study. Thus, e 2 and w 2 are not variables, but just 
parameters of the problem. Hereafter, we note them e 2 and 
■572, respectively. 

The evolution dictated by the resonant Hamiltonian, 
Eqs. (UJ) and (|6}, is obtained after applying the following 
conical transformation on the mean longitudes 



<t> = -Ai/(p + q) , 

i> = pAi — (p + q)X 2 + qm 2 , 

and on their conjugated momenta 

$ = -(p + q)Ai -pA 2 , 
* = -A 2 /(p + q) . 

In the new variables, the Hamiltonian reads 

H (V, *, hi,kx) = Ho($, *) - e^i($, *) cos V 



-9-1 



H 2 (&, ty)(ki cos(^ — zu 2 ) — hi sin('0 — ^2)) 



(8) 



(9) 



(10) 



It should be stressed that ki and hi are not conical vari- 
ables. However, in the limit of vanishing eccentricity ei, the 
equations of motions are greatly simplified 



dH 
9$ 



$ = 



dH 



* = 



dH 
9* 



fci 



1 dH 
Aldhl 



OH ■ 1 dH 

' dtp ' 1 ~ Ai dki 



(11) 
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Since the Hamiltonian (|10l) is independent of the cyclic vari- 
able cj>, its conjugated m omentum cj> is constant. In the fol- 
lowing, we assume, as in IWisdoml l| 1980ft . that the functions 
Hi(&, ^/), i — 1,2 of the interaction part are well approxi- 
mated by Hi($, where ^> = (*} is the averaged value of 



4.1 Outside of MMR 

Outside of the p : p + q MMR, we assume that the Keplerian 
Hamiltonian 'Ho($, ^) is sufficiently well approximated by 
the linear term in its Taylor series about ^ — ty. In that 
case, the resonant Hamiltonian (|10|) becomes 



(*-*) - e^i($,*)cosV 



Hi(^>, ^0 (fei cos^ — tu2) — /ii sin(-0 — ^2)) 



with, 



#,($,*) = 2Ami—aCi(a) 
m 



(12) 



(13) 



9* 



= oj^ ($, *) = pm — (p + q)n2 = +pnie . (14) 



In <|14[) . the parameter e, which cancels out at the resonance, 
is a measure of the closeness to the exact resonance. Its 
expression is 

1 + q U2 



E = 



1 - 



111 



(15) 



For infinitely small eccentricity ei, the evolution of ip 
and $ governed by the Hamiltonian ()12|) is independent of 
fci and hi. Using (|11|) . one finds 



#(t) = *(0) + 



gi(*.?) g«| 
^($,*) 2 ' 



^(0)) 



(16) 



To get the evolution of 0, we assume that &H/d& is well 
approximated by the linear expansion of &Ho/d& about ^ = 
'i', whose the expression reads 



dHo 



+ 3 



pni * - * 



p + q p + q Ai 
Integrating ([TTJl with * = (*(t)}, one gets 
nit 



#*) = 0(0) - 



p + q 



+3 



pni e| 



(17) 



(18) 



'sun/>(i) - sin(V>(0))) . 



p + gcj^*,*) 2 Ai 
From the equations of change of variables ([8]), using the 
expression of Hi($, *), Eq. JT3J, uty(<&,#), Eq. flijl. and 
removing the linear terms in (JTHJ) , one obtains 



<5Ai = -- 



6 U12 aCi(a) _ 



p mo 



sin ?/)(t) 



(19) 



The expression of Ski is obtained from the integration of 
the equations of motion (|lip of the Hamiltonian (|12|l . After 
removing the linear terms, the result is 



Ski = =f~ 



2 77i2 aC2 (a) _, 



9 -i 



p m 



- W2 



(20) 



From the expressions (|19|) and ()20|) . one can see that the 
ratio between the amplitudes of 8X1 and 8k\ is of the order 
of e^/s. Thus, for a system sufficiently close to a MMR, 
£ i$ e2, the amplitude of the TTV signal ((3| is dominated 
by 8\\. In that case, 



6 m 2 aCi (<*)_, . 
dt ~ ^+-^ e \ sm 



(21) 



pni mo 

The period of the dominant sinusoidal oscillation in the TTV 
signal P — 2?r/ Lty (<&, is related to the period of the 
transiting planet by (see Eq. ((T4j)) 



P= — Pi 

pe 



(22) 



This relation has been inverted to compute the distance e 
from each resonance in the different simulations of the sec- 
tion [3] see table [T] 

The analytical approximation of TTVs (121[) is plotted 
as gray curves for each numerical experiment in Fig. [2] The 
agreement with the numerical simulations is very good ex- 
cept for the largest distances to MMR: e ~ 1.30. Neverthe- 
less, the order of magnitude of the amplitude is still correct 
within a factor 2. 



4.2 Inside of MMR 

The reasoning is the same as in the case outside of MMR. 
The only difference is that the linear terms in the Taylor 
series of Ho(&, about \E' = $ = ^ rcs is zero. Thus, we 
expan d Ho(&, up to the quadratic order as in IWisdoml 
(1980). Furthermore, we neglect the contribution Shi in the 
TTV, as in the previous subsection. The part proportional 
to H2($, ^rcs) is thus removed from the Hamiltonian. We 



H = Ho($,* rcs ) + 



1 d 2 H 



2 a* 2 



(*-* r 



*=*re 



-e^i?i($,1' r es)cosV> 

with 
1 d 2 H 



2 9* 2 



3 p ni mi + am,2 
2 Ai am,2 



(23) 



(24) 



The Hamiltonian (|23[) is that of a pendulum. The center of 
libration depends on the sign of Hi(&, 'I'res). If it is positive, 
as in most of the cases, the libration is around ip = ifr le3 = n, 
otherwise it is around ip rcs = 0. Since the initial conditions 
used for the simulations in section [3] (Ai = A2 = VJ2 = 
270 deg) imply ?/>(0) = 0, most of the resonances have been 
crossed at the hyperbolic point. To study the TTV inside 
of a MMR, we thus change the initial conditions, and take 
A2 = 90 deg. The other parameters are those of the set 3 in 
the tabled] 

To get the evolution of cj>(t), and then <5Ai (t) , one 
should note that the oscillating term in cj> = &H/d& (|17|l 
is proportional to (>]/ — ^rcs) which is also proportional to 
ip = dH/dty . One obtains 



<Kt) = m 



nit 
p + q 



3 pni I d 2 Ho 



Ai p + q \ 9* 



(25) 
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Figure 3. Detail of the TTV amplitude through the 3 : 7 MMR. 
The initial conditions are those of the set 3 (Tab. [TJ with the 
initial mean longitude A2 = 90 deg. Each black point corresponds 
to a different initial semi-major axis 0,2 of the perturber. The thin 
solid line corresponds to the analytical formulae 1211 1 , 127t taking 
into account the amplitude of libration Aip inside the resonance. 
The amplitude of libration is computed from the initial conditions 
using the relation H(i>(0), *(0)) = H^res ± Aip, * rcs ) where U 
is the resonant Hamiltonian 1231 1. 



Removing the linear part of (|25JI . and using Eqs. 
one finds 



5Ai = 



1 



ffffl2 



-(V(*)-iM 



p mi + am2 
The corresponding TTV expression © is 



St = 



am,2 



-(lp(t) - lp rcs ) 



(26) 
(27) 

pni mi + am,2 

This expression is similar to that given bv lAgol et ail ([2005) 
for first order MMR. In (|27p . the resonant angle is librat- 
ing. So, the absolute value of (ip{t) — iptea) is bounded by the 
amplitude of libration, which is also bounded by n. 

Figure [3] shows the evolution of the dispersion a {TJ of 
the TTV signal for systems in the vicinity the 3 : 7 MMR. 
The analytical approximation plotted as a solid curve is 
given by flU} outside of the MMR, and by (J57| inside of 
the MMR, within (a 2 /a rcs ) 3/2 ~ 0.995 and 1.004. In both 
cases, the analytical values are divided by y/2 to convert 
amplitudes into dispersions. Outside of the MMR, the dis- 
persion is larger for systems closer to the separatrices. This 
is due to the factor 1/e 2 in pi"). Inside of the MMR, the 
dispersion of TTV is also larger towards th e separatrices, 
but the reason is different (jVeras et al.ll201ll ). Inside MMR, 
the dispersion is proportional to the amplitude of libration 
which is maximal at the separatrices. 

For completeness, we provide also the libration fre- 
quency ujib at the center of resonant islands. This is an 
upper limit of the frequency of the highest amplitude oscil- 
lation in TTV signals. Indeed, closer to the separatrices, the 
period gets longer. The expression of this frequency is 



wiib 



d 2 H 



d^ 2 



(28) 



; prm/6 I a— + — ) \Ci(a)\ e\ 
■m mo/ 



5 CONCLUSION 

A terrestrial planet at a period ratio larger than 2 with a 
Jovian transiting planet may produce detectable TTV only 
if the system is at resonance, but not in the exact center, 
or near period commensurability. If a system is sufficiently 
close to a second or higher order resonance, to be detected, 
then the TTV signal is dominated by a sinusoidal oscilla- 
tion which is difficult to distinguish from one resonance to 
an other. Moreover, the derived analytical approximations of 
TTV signals Eqs. (|21|) and (|27|) show that once the main pe- 
riod P of the TTV is known, the amplitude still depends on 
the mass m2 , the eccentricity e 2 , and the period P2 (through 
the semi-major axis ratio a) of the perturber. The problem is 
thus highly degenerate. The examples presented in this letter 
show that an Earth mass planet is able to produce a signal 
comparable to that produced by a Saturn mass planet. We 
thus conclude that the TTV method may benefit from radial 
velocity observations in order to characterize non-transiting 
planets. 
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